Slow fluctuations in enhanced Raman scattering and surface roughness relaxation 
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We propose an explanation for the recently measured slow fluctuations and "blinking" in the 
surface enhanced Raman scattering (SERS) spectrum of single molecules adsorbed on a silver col- 
loidal particle. We suggest that these fluctuations may be related to the dynamic relaxation of the 
surface roughness on the nanometer scale and show that there are two classes of roughness with 
qualitatively different dynamics. The predictions agree with measurements of surface roughness 
relaxation. Using a theoretical model for the kinetics of surface roughness relaxation in the pres- 
ence of charges and optical electrical fields, we predict that the high-frequency electromagnetic field 
increases both the effective surface tension and the surface diffusion constant and thus accelerates 
the surface smoothing kinetics and time scale of the Raman fluctuations in manner that is linear 
with the laser power intensity, while the addition of salt retards the surface relaxation kinetics and 
increases the time scale of the fluctuations. These predictions are in qualitative agreement with the 
Raman experiments. 



In a recent experiment by Weiss and Haran [1], large 
spectral fluctuations in the relative intensities of differ- 
ent Raman lines that varied on a time scale of a few tens 
of seconds, were measured in the surface-enhanced Ra- 
man scattering (SERS) of single rhodamine 6G molecules 
adsorbed on silver nanocrystals. The rate of spectral 
fluctuations was demonstrated to increase with laser in- 
tensity and decrease with addition of salt in the solution 
of the silver colloidal particles. In addition, a decay of 
the overall intensity of the scattering was observed on 
a scale of the order of hundreds of seconds; this decay 
was also correlated with the laser intensity [1]. Finally, 
fluctuations of the overall intensity on a time scale com- 
parable with that of the fluctuations of the individual 
spectral lines were measured in Refs. [I]. These fluctu- 
ations, sometimes termed "blinking", are frequently ob- 
served in single-molecule SERS studies [2-5]. The long 
time scale observed in the modulations of the spectrum 
of molecules adsorbed on silver colloidal particles, was 
suggested [1] to arise from (slow) motion of the adsorbed 
molecule, which leads to variation of a charge-transfer 
interaction between the molecule and the surface. 

In this paper, we propose an alternative interpreta- 
tion of both the slow SERS spectrum fluctuations, and 
of the decay and the fluctuations of the overall intensity 
("blinking") of the SERS spectrum measured in Ref. [1] 
and show that the time scales for these phenomena are 
consistent with the relaxation of the surface roughness 
on the nanometer scale [6] . It is well known that a signif- 
icant part of the enhancement of the SERS signal is due 
to the surface roughness [7,8]. Using theoretical mod- 
els for the kinetics of surface roughness relaxation in the 



presence of charges and optical electrical fields, we pre- 
dict the dependence of the time scale of the fluctuations 
on the laser power intensity and on the amount of added 
salt in agreement with the observations. 

Conventionally, the SERS enhancement mechanism is 
separated into the electromagnetic (EM) field enhance- 
ment and the chemical enhancement [3] . The relaxation 
of surface roughness can affect both the EM and chemi- 
cal mechanism of the SERS enhancement, depending on 
the characteristic length-scale and the amplitude of the 
relaxation modes. The surface relaxation dynamics it- 
self, however, is governed by the EM mechanism in our 
model, as we elaborate below. These constitute the start- 
ing point of our analysis. 

Our dynamical model applies both to equilibrium fluc- 
tuations and to the decay to equilibrium (a smooth sur- 
face) of non-equilibrium, surface roughness that can have 
large amplitudes, depending on the initial preparation of 
the nanoparticles. The amplitude of equilibrium fluc- 
tuations is typically small for perfectly smooth surfaces 
[9]. However, rough surfaces (where the roughness due 
to sample preparation decays slowly) show larger thermal 
fluctuations [10]. Thus, the roughness fluctuations may 
have a strong effect on the SERS enhancement mecha- 
nism and may lead to the observed fluctuations in the 
relative intensities of the different Raman lines. 

We first estimate the time scale for relaxation of the 
surface roughness. We use a model for surface relaxation 
developed by Mullins [11] based on an isotropic expres- 
sion for the surface energy. This theory is applicable 
either above the roughening transition temperature or 
for vicinal surfaces (that are not perfectly smooth on the 
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atomic scale), even below the roughening temperature. 
The dynamics of ideal, high symmetry surfaces involves 
the creation of steps and facets [12] and is outside the 
scope of our work. Assuming that surface diffusion (with 
diffusion constant D s for single atom motion on the sur- 
face) is the only relevant process, leads one to the kinetic 
equation: 



where h(p,t) is the local height, 70 is the surface ten- 
sion, p — (x,y) is the in-plane position vector and V = 
(d x , d y ), and T = D s a 4 jo / kgT with a being an atomic 
length scale (e.g. the nearest-neighbor distance, for sil- 
ver a ~ 2.8 A). To estimate the characteristic time of the 
SERS spectral fluctuations, we analyze the height-height 
dynamic correlation function, Chh(p, t) = (h(p, t)h(0, 0)) . 

Chh(p,t) can be straightforwardly obtained from Eq. 
(1) (see e.g., Refs. [9,13]) assuming that the surface 
atoms are in contact with a thermal bath: Chh(p,t) = 

few -C ^-J—MQP) d <7, where J (x) is the Bessel func- 
tion, and q = 2tt/L 7 with L being the system size. To 
estimate the characteristic relaxation time, r , numeri- 
cally, one can adopt the criterion: Chh(po, i~)/Chh(po, 0) = 
1/ exp(l), were po is of the order of a size of the adsorbed 
molecule. 

An accurate estimate of the relaxation time, r, requires 
a knowledge of the surface diffusion constant, D s , of sil- 
ver estimated in the range D s ~ 1.8 x 10~ 15 cm 2 / sec 
[14] to D s ~ KT 14 cm 2 / sec [15]. Using L ~ 50 nm of 
the order of a size of silver colloids in the experiment [1] , 
surface tension of Ag [16], 70 — 1500 erg/cm 2 , po ~ 6 
A, we find that the correlation time of the equilibrium 
fluctuations is r ~ 18 sec. This shows that the time 
scale for the equilibrium relaxation of surface roughness 
on nanometer scale at room temperature is comparable 
to the correlation times of the slow fluctuations of the 
SERS spectrum measured in [1] . An analytic estimate for 
the correlation time can be obtained from the asymptotic 
long-time form: C hh (p ,t) ~ C hh (0,t) ~ ffi ""fc^ , 
where r = 1/Tq 4 . We emphasize that the relaxation 
time, r , is sensitive to the numerical values used for the 
lattice constant as well as to the wavelength associated 
with the size of a colloid. 

In addition to the small amplitude, equilibrium fluc- 
tuations of the surface, much larger amplitude surface 
roughness may arise due to the sample preparation. 
These features tend to decay in time, leading to a smooth 
surface at the atomic level and this relaxation to equilib- 
rium of a relatively smooth surface (and hence less en- 
hancement of the Raman scattering), leads to the over- 
all decay of the SERS spectrum intensity. The loss of 
the Raman signal is a well known phenomenon that also 
occurs in electrochemical [17,18] and ultrahigh vacuum 
[19,20] SERS systems. This effect was attributed to a 



diffusive loss of surface adatoms in a number of investi- 
gations (see e.g., introduction in Refs. [17,18]); however, 
none of these, to the best of our knowledge, explain the 
slow time scale in terms of cooperative surface tension 
effects nor do they discuss the dependence of the decay 
on the laser optical field and on added salt. 

Another type of experiment to which our kinetic model 
is relevant, directly measures the dynamics of artificially 
created, nanoscale surface roughness features [15]. Since 
the surface profiles in the experimental systems are, in 
general, far from being simple sinusoids, we predict the 
relaxation time for two important classes of surface fea- 
tures. 

To estimate the relaxation time for the two classes of 
surface features we consider two specific surface rough- 
ness profiles at an initial time t = 0: (i) a non-mass- 
conserving profile, modelled as a Gaussian protrusion 
where ft > 0, h(p, t = 0) = ho e~ ap and (ii) a mass- 
conserving profile, modelled as a region where a protru- 
sion with ft > is adjacent to an indentation where 
ft < 0. The average value of ft over the entire surface 
is zero and for convenience we consider: h(p, t = 0) = 
fto e~ Qp (1 — ap 2 ), where or 1 / 2 is the characteristic lat- 
eral scale of the feature, and ho is the amplitude. The ex- 
periments of Ref. [15] correspond to the mass-conserving 
case (ii) (possibly because of the manner in which the 
surface was scratched) , while we expect that general sur- 
face roughness of colloidal particles, applicable to the 
Raman experiments, to be more similar to the non-mass- 
conserving case (i) . It is straightforward to compute the 
time evolution of these two types of profiles, by solving 
Eq. (1) with the corresponding initial conditions of cases 

(i) and (ii). The time dependence of the decay of the 
maximum height (located at the origin, p = 0) of these 
Gaussian-like peaks can be obtained analytically. For 
case (i), we find h(0,t) = h (nTg/t) 1 / 2 (3(T g /t) . For case 

(ii) we find ft(0, t) = ho [2 T g /t - 2^( Tg /t) 3 / 2 0( Tg /t) ] , 
where (3(x) =cxp(x)crfc(y / a;) and T g — ksT/ (64-f D s a 2 a 4 ) 
is the characteristic decay time, and erfc(x) is the com- 
plementary error function. 

The most important observation is that these two types 
of profiles have qualitatively different smoothing kinet- 
ics - the mass conserving profile (ii) decays much faster 
since in this case the transport of matter need only oc- 
cur near the boundary between the protrusion and the 
indentation; that is, atoms are locally transferred from 
the region where ft > to the region where ft < 0. 
For case (i) of the non-mass-conserving profile, the mat- 
ter must be transported to a much larger scale. This 
is reflected in the expressions for the asymptotic:, long 
time (r g <C t) evolution of the height. For case (i), 
the decay of the maximum of the profile (located at 
the origin), h(0,t) ~ ft n 1 / 2 (Tg/t) 1 / 2 . For case (ii) 
h(Q,t) ~ 2 ft T g /t tends to zero much faster. 

We now use this model for surface smoothing in case 
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(ii) to estimate the decay time for the (approximately) 
mass-conserving surface features studied in Ref. [15]. Us- 
ing their estimate for the diffusion constant at zero ex- 
ternal potential, D s ~ 10~ 14 cm 2 / sec and a value for 
0,-1/2 ^ g nm tnat correS p 0n d s to the extent of the 

scratch, we obtain that a mass conserving, Gaussian pro- 
file with initial amplitude ho ~ 2 nm decays to an atomic- 
scale estimated as h(r) ~ 0.3 nm, in a time r ~ 200 sec. 
This theoretical estimate is consistent with the exper- 
imentally measured times in the scratch experiments of 
Hirai [15]. An estimate of the decay time for case (i) (i.e., 
the non-mass-conserving profile) , yields for the same pa- 
rameters, a decay time that is about 15 times slower, 
t ~ 3100 sec. In the case of colloidal particles, we ex- 
pect that typical, non-equilibrium features (preparation 
dependent roughness) may be non-mass-conserving pro- 
trusions or indentations and that the kinetics of case (i) 
would apply. We estimate, for instance, that a more lo- 
calized, Gaussian protrusion of amplitude 1 nm and ex- 
tent or 1 / 2 ~ 6nm would decay to an atomic size of 0.3 
nm in a time 200 sec, consistent with the times measured 
in the Raman experiments [1]. 

The Raman experiments show a systematic depen- 
dence of the modulation and relaxation time scales on 
the salt and the electric field as mentioned at the begin- 
ning. Since the characteristic time scales vary inversely 
with the product of the surface tension and diffusion con- 
stant (see Eq. (1)), we consider: (i) the effect of salt on 
the surface tension (ii) the effect of the laser field on the 
time scale via its effect on the surface tension (iii) the 
effect of the laser field on the surface diffusion constant. 
To treat the effects of surface charges [21] and the laser 
electric field on the silver surface tension and hence on 
the surface relaxation, we model the charged colloidal in- 
terface as an elastic, almost planar surface, with a fixed 
and uniform surface density of charge a, and height h(p). 
The surface tension, given by calculating the free energy 
cost (including the electrostatic effects due to the charges 
and the salt) of deviations of the surface from the planar 
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the valence of salt ions, and n is the salt concentration, 
E s = 4ir<r/e is the electric field at the colloid surface. 
This result can also be obtained from a scaling argument: 
the tension is the product of the energy density and the 
volume divided by the cross-sectional area; this is pro- 
portional to energy density eE 2 multiplied by a charac- 
teristic length, which here is the Debye length, 

This result shows that with the addition of salt, the ef- 
fective surface tension decreases, 7 e ~ nT x l 2 . This leads 
to an increase of the relaxation time, r, as observed. Tak- 
ing the experimental values of the parameters used above, 
we see that in order to obtain an effect of order unity (i.e., 
an effective electrostatic surface tension, 7e , equal to the 



bare surface tension of silver, 7e — 70 ~ 1500 erg/cm 2 ) 
the surface of a silver colloidal particle must have a sur- 
face charge density, a ~ 10e/nm 2 in 100 mM salt solu- 
tion. This is a rather high surface charge density, but is 
still within a realistic range for the experiments. 

To quantify the contribution of the laser field to the 
surface energy and diffusion constant of a silver nanocrys- 
tal requires an accurate value for the enhanced surface 
electric field; this requires an accurate model of the mi- 
croscopic mechanism of the enhancement, which is not 
yet completely understood [7]. However, we can predict 
the functional dependence of the both the surface ten- 
sion as well as the surface diffusion constant on the field 
at the surface. The surface energy density, u, of a semi 
infinite metal sample with a planar surface boundary in 
the presence of a high-frequency, optical electromagnetic 
field, is [23]: 



u(z) 



1 

8^ 



d(cjs(u))) 

diO 



\E {z) \ 2 + d ±^\H {z) \ 2 



(2) 



where E(z) = E e tkz is the spatially-dependent part 
of the electric field vector with amplitude, Eq, and the 
plane, z = 0, corresponds to the metal interface plane; 
the same definitions apply to the magnetic field, H. 
For the frequency-dependent dielectric function, e(u), we 
use the simplest dispersion model of the free, classical 
electron gas [23]: e(u>) = 1 — (^f) 2 , where u p ~ 9.2 
eV~ 1.4 x 10 16 sec -1 is the plasma frequency of silver, 
and u> ~ 3.5 x 10 15 sec -1 for the 530 nm wavelength 
laser beam. For the optical frequencies relevant to the 
experiments [1], fi{un) = 1. Using the dispersion relation 



and 



[23], k 2 = ^w 2 , we obtain: k 
therefore, both the electric and magnetic fields decay ex- 
ponentially within the silver colloidal particle with a typ- 
ical decay length [24], x _1 , of 22 nm. In this frequency 
region therefore, the metal is reflecting with a skin depth, 
X -1 . Using the relationship between H and E [23], we 
obtain the time-averaged energy per unit surface area 
{i.e., effective surface tension), -f em = J Q °° u(z) dz, in the 



presence of an electromagnetic field: j e 



8-n X 



Again, 



this can be predicted by a scaling argument: the ten- 
sion is the product of the energy density (proportional 
to the laser intensity) and the characteristic length, here 
the skin depth, in qualitative agreement with the exper- 
iment. 

However, an estimate of the magnitude of this effect us- 
ing the bare values E ~ 0.6(erg 1 / 2 /cm 3 / 2 ) ~ 180 V/cm, 
(corresponding to a laser power density of 100 W/cm 2 ) 
and x^ 1 ^ 22 nm is 7em ~ 3 x lO^ 8 ^. This quan- 
tity, 7 em , is thus about 10 11 times smaller than the bare 
coefficient of surface tension of silver, 70 ! It is known 
that the electric field at the surface in the case of SERS 
is greatly enhanced by the surface plasmon resonances in 
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the system [7,25]. Recent experiments [2,3] report that 
the SERS enhancement factor, (theoretically predicted to 
be proportional to (£/E ) 4 ) where £ is the intensity of 
scattered, enhanced optical field (see also Refs. [27,28])), 
reaches a value of 10 14 — 10 15 . This enhancement is still 
not large enough to bridge the gap between our predic- 
tion for the optical field-induced effective surface tension, 
7em, and the bare surface tension of silver, 70. We note, 
however, that the intensity of the localized field at the 
surface of a metal colloid may exceed the intensity of 
the radiated, enhanced optical field [7]. Atomic scale, 
globally distributed roughness (adatoms, terraces, kinks, 
small islands) may produce an extra enhancement of the 
electric field at the surface of a colloidal particle, as well, 
although a quantitaive understanding of this mechanism 
is still lacking (see e.g., [26]). In summary, the scaling 
of the tension and hence the time scale for surface re- 
laxation with the electric field is qualitatively correct; 
however, the magnitude depends on the actual value of 
the surface field. 

Although the effect of the field on the surface tension 
may turn out to be too small to matter, the change in 
the local diffusion constant of silver adatoms with field 
may be more significant; the diffusion constant enters the 
kinetic equations discussed above. Indeed, in the case of 
a static electric field applied to a roughened metal inter- 
face, Hirai et al. [15] show from their experiment that 
the surface diffusion coefficient, D s , depends exponen- 
tially on the applied potential difference. This is in ac- 
cord with theoretical calculations [29] of surface diffusion 
in the presence of external fields, in the limit where the 
electrostatic energy is much larger than the thermal en- 
ergy. This is the case in the static experiments where 
the potential drop occurs on a length scale of ~ 1 nm of 
order of the Debyc length; this gives rise to a very large 
electric field. However, this limit is not applicable to the 
Raman experiments. The electrolyte cannot respond to 
the high frequency optical field and provides no screen- 
ing of the laser field. The potential drops over a length 
scale given by the optical wavelength, and not the Debye 
length. The resulting field is thus orders of magnitude 
smaller than in the static case. Indeed, the Raman ex- 
periments show only a linear dependence of the rate of 
the SERS spectral fluctuations on the laser intensity and 
no exponential behavior. 

Calculations [29] of surface diffusion in the presence 
of external fields in the limit where the field energy is 
smaller than the thermal energy, show that the surface 
diffusion coefficient depends linearly on the laser inten- 
sity [30]: D s ~ D S {E = 0) (l + W 2 ) where W = 
eE a/(4:kBT) <C 1. This trend agrees with the exper- 
imental observation that the time scale varies with the 
laser intensity. However, as in the case of the surface ten- 
sion, this requires a significant enhancement of the local, 
surface field to yield a measurcablc effect. We note that 
even if the field contribution to the tension is small, the 



effect on the diffusion constant may be significantly larger 
since the field dependence of the tension depends on the 
ratio of the intensity to the bare silver surface tension 
while the correction to the diffusion constant varies with 
the ratio of the intensity to the thermal energy, which is 
about two orders of magnitude smaller than the energy 
associated with the surface tension of silver. 

The best way to verify our predictions would be to per- 
form SERS measurements on surfaces with in situ con- 
trol of surface roughness by e.g., the methods described 
in Rcf. [15]. Monitoring the SERS spectrum as a func- 
tion of the surface roughness, the salt concentration, and 
the laser intensity should provide an ultimate test of our 
idea. 
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